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ABSTRACT 

The presence of fine structures in the sunspot vector magnetic fields has been 
confirmed from Hinode as well as other earlier observations. We studied 43 
sunspots based on the data sets taken from ASP/DLSP, Hinode (SOT/SP) and 
SVM (USO). In this Letter, (i) We introduce the concept of signed shear angle 
(SSA) for sunspots and establish its importance for non force-free fields, (ii) We 
find that the sign of global a (force-free parameter) is well correlated with the 
global SSA and the photospheric chirality of sunspots. (iii) Local a patches of 
opposite signs are present in the umbra of each sunspot. The amplitude of the 
spatial variation of local a in the umbra is typically of the order of the global a 
of the sunspot. (iv) We find that the local a is distributed as alternately positive 
and negative filaments in the penumbra. The amplitude of azimuthal variation of 
the local a in the penumbra is approximately an order of magnitude larger than 
that in the umbra. The contributions of the local positive and negative currents 
and a in the penumbra cancel each other giving almost no contribution for their 
global values for whole sunspot. (v) Arc-like structures (partial rings) with a 
sign opposite to that of the dominant sign of a of the umbral region are seen at 
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the umbral-penumbral boundaries of some sunspots. (vi) Most of the sunspots 
studied, belong to the minimum epoch of the 23 rd solar cycle and do not follow 
the so-called hemispheric helicity rule. 

Subject headings: Sun: magnetic fields, Sun: photosphere, Sun: sunspots 



INTRODUCTION 



Helical patterns in the sunspots and associated features have been observed for a long 
time (jHald Il925l . 119271 ; iRichardsonl Il94ll ) with a hemispheric preference of their chirality, 
which is independent of solar cycle. Since the 90's, the subject has been intensively revis- 
ited and the similar behaviour of hemispheric patterns fo r various solar f eatures have been 



report e d by many researcher s (IHagino fc Sakurail (12004 ): iNandyl (120061 ): iBernasconi et al. 



( 120051 ): iPevtsov fc Longcopd (120011 . 120071 ) and references therein). However, this hemi- 
spheric behaviour needs further investigation due to some inconsistencies reported for dif- 
ferent phases of a solar cycle and also for da ta sets obtained from different magnetographs 
( IHagino fc Sakuraill2005l : IPevtsov et al.ll2008h . 



For a force-free field, the glo bal twist per unit a xial length is given by the force-free 



parameter a (see Appendix A of iTiwari et al.l (120091 )). Some recent studies ( ITiwari et al. 



20081 . l2009al fbh have shown that the global a of an active region bears the same sign as its 
associated features/structures observed at chromospheric and coronal heights. The chromo- 
spheric and coronal sign of twist is inferred from the topological chirality sign of the observed 
features. This leads us to believe that some form of the photospheric global twist exists on 
the scale of sunspots. However, the structures in the sunspot fields revealed by modern 
vector magnetographs with high spatial and spectral resolution compels us to make a careful 
revaluation of global a and its physical meaning. 



Since the photospheric field is not force-free (IMetcalf et al.lll995l ). we need an alternative 
measure of the twist other than a. We introduce the concept of signed shear angle (SSA) 
for sunspot magnetic fields in this paper and show how SSA is directly related to vertical 
current (J z ) and a, irrespective of the force-free nature of the sunspot fields. 

The pre s ence o f oppositely directed currents in a single unipolar sunspot was first shown 
by ISevernyil (119651 ) . For a detaile d investigation of lo cal a distri bution in three sunspots 
using 46 vector magnetograms see IPevtsov et al.l (Il994h . Recently, ISu et al.l (120091 ) reported 
an interesting pattern of fine structures in the a distribution within one active region (AR) 
using Hinode data with higher resolution. We present a comprehensive study of 43 sunspots 
with high resolution and establish the contribution of such fine structures to the global twist. 
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For this purpose we will rely on J z and a values. 



The helicity hemispheric rul e or more precisely twist hemispheric rule is claimed to be 



established by many researchers (Seehafer 199ol; Pevtsov et 



Bao & Zhan 5ll998[ lLongcope et al.lll998l: lHagino & Sakuraill2005h and has recently been 



1995; Abramenko et al. 1996 



matter of some d e bate (jHagino fc Sakuraill2005l ; iPevtsov et al.ll2008l ). A model developed by 
Choudhuri et al.l (120041 ) predicts deviation from the twist hemispheric rule in the beginning 
of the solar cycle. However, some observers claim this devi ation from the hemis pheric rule 
may be present in different phases of different solar cycles (jPevtsov et al.ll2008l ). We have 
studied forty three AR's (as shown in Table 1) mostly observed during the declining phase 
of the solar cycle 23. None but five follow the twist hemispheric rule. 



In the following section (Section 2), we discuss the data sets used. Section 3 describes 
the analysis and results obtained. Finally in Section 4 we present our conclusions. 



DATA SETS USED 



We have used the vector magnetograms obtained from So l ar Op t ical Telescope /Spectre - 



polarimeter fSO T/SP: iTsuneta et al.l (120081 ) ; ISuematsu et all ( I2008f ) ; llchimoto et al.l (120081)') 
onboa rd Hinode ( Kosugi et al. 2007 ) and the Advanced Stokes Polari meter (ASP;lElmore et al. 



( 1992 ) ) as w ell as the Diffraction Limited Spectro-polarimeter (DLSP; ISankarasubramanian et al. 



( 120041 . 120061 )) of the Dunn Solar Telescope (DST). Standard and well-established calibration 
procedure was adopted for ASP/DLSP data. Th e procedure for obtaining the vector fields 
from the A SP/DLSP data are described elsewhere (jElmore et al. Il992l ; ISankarasubramanian et al. 
2004J . l2006h . 



The Hinode (SOT/SP) data have been calibrated by the standard "sp_prep" routine 
available in the Solar-Soft packages. The prepared polarization spect ra have been inverted to 
obtai n vector magnetic field components using an Unno-Rachkowsky ( UnnolE956 : Rachkowskv 



1982 



1967 ) inversion under the assum ption of Milne-Eddington (ME) atmosphere (ILandolfi fc Landi Degl'Innoce: 



Skumanich fc Liteslll987). The 180° azimuthal ambiguity in the data sets are rem oved 



by using acute angle method (lHarveylll969l ; ISakurai et al.lll985t ICuperman et al.lll992l ). All 
the data sets used, have high spatial sampling. For example ASP ~ 0.3 arcsec/pix, DLSP 
~ 0.1 arcsec/pixel and Hinode (SOT/SP) ~ 0.3 arcsec/pixel. However, a few observations 
are seeing limited to about an arcsec. 

To minimize noise, pixels having transverse (B t ) and longitudinal magnetic field (B z ) 
greater than a certain level are only analyzed. A quiet Sun region is selected for each sunspot 
and la deviation in the three vector field components B x , B y and B z are evaluated separately. 
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The resultant deviations in B x and B y is then taken as the la noise level for transverse field 
components. Only those pixels where longitudinal and transverse fields are simultaneously 
greater than twice the above mentioned noise levels are analyzed. 

The data sets with their observation details are given in Table 1. The data sets observed 
from August 2001 to April 2005 are obtained with the ASP and those observed from June 
2005 to December 2005 are from the DLSP. Two vector magnetograms observed on 09 
January 2007 and 06 Fe bruary 2007 f r om Solar V ector Magnetogra ph at Udaipu r Sola r 
Observatory (SVM-USO: iGosain et all (booi l2006h ) and reported in iTiwari etaD J2008h . 
have also been included to improve the statistics. All the other data sets obtained since 
November 2006 onwards are taken from Hinode (SOT/SP). 



DATA ANALYSIS AND RESULTS 



We have used the following formula to compute the local a values 

(V x B) z 



a 



B-, 



(1) 



T he global a va l ue of the active regions is estimated from the following formula as described 



m 



Tiwari et all (120091 ) 



On 



9By 

> dx 



(2) 



This estimate was shown to be not seriously affected by the polarimetric noise ( Tiwari et al 



20091 ). Moreover, since a g is weighted by strong field v alues (|Hagino fc Sakurai 



20041 ) and 



not affected by singularities at polarity inversion lines (ITiwari et al.l 120091 ) . this parameter 
is more accurate than a simple average of local a. 



Hagyard et al.l (119841 ) introduced the shear angle A$ = $ 



obs 



$ p0 i, where § obs and 



$ pot are the azimuthal angles of observed and potential fields respectively. The amplitude of 
this angle was studied at the polarity inversion lines to investiga te the flare related changes 
(IHagyard et al.lll990l ; lAmbastha et al.lll993l ; lHagyard et al.lll999l ). To emphasize the sign of 
shear angle we wish to introduce the signed shear angle (SSA) for the sunspots as follows : 
We choose an initial reference azimuth for a current-free field (obtained from observed line of 
sight field). Then we move to the observed field azimuth from the reference azimuth through 
an acute angle. If this rotation is counter-clockwise, then we assign a positive sign for the 
SSA. A negative sign is given for clockwise rotation. This sign convention will be consistent 
with the sense of azimuthal field produced by a vertical current. This sign convention is also 
consistent with the sense of chirality (for details, see Appendix A). The potential field has 
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been computed using the method of ISakurail (119891 ) . Mean of the SSA obtained for a whole 



sunspot is taken as the global value of SSA for that sunspot. 

The force-free parameter a involves three dimensions since it basically represents the 
rate of change of rotation per unit axial length. The SSA is the rotational deviation of the 
projection of the field onto the photosphere from that of a reference current-free field. The 
a parameter is a gradient of angle per unit length while SSA is just an angle. We therefore 
cannot expect a strong correlation between the amplitudes of both the quantities, the SSA 
and the a parameter. But we do find a good correlation between their signs as evident from 
Table 1. 

Signed shear angle (SSA) provides the sign of twist irrespective of whether the photo- 
spheric magnetic field is force-free or not. Table 1 shows that the sign of a is same as the 
sign of SSA. Thus, we conclude that even if the photosphere is non force-free, the sign of 
global a will empirically give the sign of global SSA and therefore the sign of global twist 
(chirality) of the sunspots. 

To avoid any kind of projec tion effect we have transformed the data sets to disk center 



(jVenkatakrishnan Sz Gary! Il989l ) if the observed sunspot is equal to more than 10 degrees 
away from the disk center. In some active regions both the polarities are compact enough 
to be studied separately. We have treated each pole of those active regions as an individual 
sunspot and this is denoted in the Table 1 after the NOAA no. of sunspots by plus or minus 
sign. 

Two examples of the local a distribution for the data sets obtained from Hinode 
(SOT/SP) are shown in Figure 1. The positive/negative contours are shown in red/blue 
colors. The local a patches are seen in the umbra and filamentary distribution of a is ob- 
served in the penumbral region. We find that the inclination angles oscillate between ~ 
30-80 degrees when we go along azimuthal direction in the filamentary pe numbral struc- 



tures . This is consistent with the interlocking-comb penumbral structure (jlchimoto et al. 



20071 ) of the penumbral magnetic fields. The vertical current J 2 has two components, viz. 



— r^T an d r~^7^' ^ we a PP roxrma te the observed transverse field (B t ) to be mostly 
radial (B t ~ B r ) then we can interpret the azimuthal variation of J z to result from the term 
_idB^ This term is not expected to contribute to global twist. However l 9 ^^ could be 

r oq> 1 ° r or 

an important contributor to the global twist. A detailed investigation of this interesting pos- 
sibility is deferred to another paper. For the present, we obtain positive and negative values 
of current side by side in the penumbra. Because the a parameter depends on the current, 
this oscillation in the filamentary structure across the penumbral filaments is expected for 
the a values too. 
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The distribution of vertical current and local a in the penumbra show higher values 
than that in the umbral regions. An arc and a straight line, selected respectively in the 
penumbra and umbra of AR NOAA 10933 have been over plotted as shown in the left panel 
of Figure 1. The corresponding values of vertical current and a along the arc and the line 
are shown in the Figure 2(a) and Figure 2(b) respectively. We can see that both the positive 
and negative vertical current as well as a are equally distributed in the penumbra along the 
azimuthal direction. This gives a negligible contribution to the global current and global a 
values thereby indicating that the contribution of is indeed small. We have selected 

an arc rather than the complete circle because many times sunspots are not circular and 
therefore selecting a proper penumbral region is not possible by a full circle. Similar arcs 
have been selected in the other sunspots and all the time it is seen that both the positive 
and the negative vertical current as well as a are distributed equally in the penumbra giving 
negligible contribution to their global values. While current and a variations are correlated 
for positive B 2 , they will be anti correlated for negative B z . 

Figure 2(b) shows a typical profile of spatial variations of current and a across the 
umbra (along the line) in the AR NOAA 10933 shown in the left panel of Figure 1. We see 
that the amplitude of variation of a in the umbra is smaller than that in the penumbra by 
approximately an order of the magnitude and is of the same order as that of the global a 
value of the whole sunspot. The current values are not as much smaller than that in the 
penumbra but are less balanced out than in the penumbra. 

In the right panel of Figure 1 an arc like structure (partial ring) with bunch of red 
contours (positive a) can be observed. This is opposite to that of the dominant negative 
global a of the sunspot. Such partial rings with opposite sign of the global value are observed 
in 10 of the sunspots from our sample. In the rest of the sunspots mixed current and a are 
present in the umbra with one dominant sign and no such specific structures are seen at the 
umbral-penumbral boundaries. 

A few sunspots in our data sets studied are small and have no penumbra. Some ASP 
data doesn't show fine structures in penumbra due to lack of spatial resolution. We have 
included these active regions in our study to look at the hemispheric behaviour of the global 
twist. 

Most of the data sets we studied are observed during the declining minimum phase of the 
solar cycle 23. All except five of the sunspots observed do not follow the twist hemispheric 
rule. 
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4. DISCUSSION AND CONCLUSIONS 

We have introduced the concept of SSA for sunspots and further find that the SSA and 
the global a value of sunspots have the same sign. Thus, a gives the same sign as the SSA 
and therefore the same sign of the photospheric chirality of the sunspots, irrespective of the 
force-freeness of the sunspot fields. As can be observed from the Table 1, the magnitudes 
of SSA and a g are not well correlated. This lack of correlation could be due to a variety 
of reasons: (a) departure from the force-free nature (b) even for the force-free fields, a is 
the gradient of twist variation whereas SSA is purely an angle. The missing link is the 
scale length of variation of twist. The magnitude of global SSA, therefore, holds promise 
for characterizing global twist of the sunspot magnetic fields, irrespective of its force-free 
nature. 

Patches of vertical current and a with opposite signs are present in the umbra of each 
sunspot studied. Since opposite currents r epel, the exist ence of a dominant current may be 



an useful binding force for the umbra (cf. iParkerl (119791 )). This will be examined in detail 



in further studies of evolution of twist in decaying sunspots. One sign of a dominates in 
the umbra which is also seen to be the sign of the global a of the sunspot. Similarly, the 
magnitude of the global a is of the same order as the amplitude of the local a in the umbra. 

The filamentary structures of vertical current and local a are observed in the penumbra 
of the sunspots and are, as discussed above, due to oscillatory behaviour in the inclination 
values and therefore gradients of the transverse field in the azimuthal direction. We find that 
the contributions of both positive and negative values of vertical current and a to their global 
values cancel each other in the penumbra of the sunspot. Thus the penumbral fine structures 
provide a negligible contribution to the global a and current values of sunspots. The mutual 
repulsion of opposite currents also seem to balance out in the azimuthal direction. It is to 
be seen whether disruption of this balance leads to sunspot rotation and change in global 
twist. At any rate, the observed balancing of the filamentary currents in the azimuthal 
direction may be an important contribution to the force-free nature of the sunspot fields. 
The amplitude of a variation is approximately an order of magnitude smaller in the umbral 
regions than that in the penumbral regions and is of the order of the global a of the whole 
sunspot. 

Partial rings with opposite signs to that of dominant sign of umbral a are observed at 
sunspot umbral-penumbral boundaries in 10 out of 43 sunspots studied. However even in 
these few cases, the rings are never complete in any sunspot. 

Most of the AR's observed do not follow the twist hemispheric rule. This issue of 
hemispheric rule needs to be reinvestigated over a longer period as well as with improved 
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data. 



Some researchers have tried to pred i ct fla re activity from loca l distribution of a as well 
as global a v alues of sunspots (INandyl (120081 ); lHahn et al.l (120051 ) and references therein). 
Nandyl (120081 ) concluded from a study of AR 6982 that the global twist present in a sunspot 
does not influence the flaring activity. It is rather, governed by the spatial distribution and 
evolution of twisted sub-structures present in the sunspot. This conclusion indeed needs 
more study. We plan to address, in our forthcoming study, the question of relation between 
flaring activity and the role of global as well as the local twist present in a large number of 
sunspots. 

For the present, we demonstrate that the sign of the SSA provides the sign of the 
photospheric chirality irrespective of its force-free nature. The sign of the global a of a 
sunspot is determined by the dominant sign of umbral a values without much contribution 
from the penumbral a values. 
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APPENDIX 



A. Relation between the sign of SSA and the sense of chirality 

The definition of the signed shear angle (SSA) is introduced in the Section 3. Figure 3 
shows four structures, first two having positive chirality and the next two having negative 
chirality. The sign of B pof and B & s point inward for negative B 2 and outward for positive 
B 2 . The rotation from B pot to B b. s through an acute angle is counter clockwise for the cases 
of positive chirality and clockwise for negative chirality. This is consistent with positive and 
negative SSA respectively, by definition. Thus, the sign of SSA will bear the same sign of 
the chirality. 
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Table 1: List of the active regions studied. The global a value, the signed shear angle (SSA) 
and other details of the sunspots are given: 



AR No. 


Date of 


Global Alpha 


Shear Angle 


Position 


Hemispheric 


(NOAA) 


Observation 


(a 9 :/meter) 


(SSA: deg) 




Hclicity Rule 


10972 


07 Oct 2007 


-2.331 x 10" 8 


-1.085 


S05W20(t) 


No 


10971 


29 Sep 2007 


3.053 x 10~ 8 


3.214 


N03W07 


No 


10970 


05 Sep 2007 


-2.001 x 10~ 8 


-0.308 


S07W58(t) 


No 


10969 


29 Aug 2007 


-3.424 x 10~ 8 


-4.488 


S05W33(t) 


No 


10966 


09 Aug 2007 


-2.539 x 10~ 8 


-3.595 


S06E07 


No 


10963(-) 


12 Jul 2007 


-2.459 x 10~ 8 


-4.636 


S06E14(t) 


No 


10963(+) 


12 Jul 2007 


-3.440 x 10~ 8 


-4.495 


S06E14(t) 


No 


10961 


02 Jul 2007 


-5.119 x 10~ 8 


-4.973 


S10W16(t) 


No 


10960 


07 Jun 2007 


3.027 x 10~ 8 


4.486 


S07W03 


Yes 


10956(-) 


18 May 2007 


9.642 x 10~ 8 


11.595 


N02E07 


No 


10956(+) 


18 May 2007 


6.458 x 10~ 8 


5.352 


N02E07 


No 


10955 


13 May 2007 


-6.737 x 10~ 8 


-1.887 


S09W35(t) 


No 


10953 


29 Apr 2007 


-6.673 x 10~ 9 


-3.071 


S10E22(t) 


No 


10944 


03 Mar 2007 


-2.084 x 10~ 8 


-4.635 


S05W30(t) 


No 


10941 


06 Feb 2007 


-2.745 x 10~ 8 


-3.069 


S07W36(t) 


No 


10940 


01 Feb 2007 


-1.948 x 10~ 8 


-4.726 


S04W05 


No 


10939(-) 


23 Jan 2007 


-3.033 x 10~ 8 


-5.105 


S04W57(t) 


No 


10939(+) 


23 Jan 2007 


-8.289 x 10~ 9 


-0.869 


S04W57(t) 


No 


10935 


09 Jan 2007 


-2.412 x 10~ 8 
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pixels pixels 

Fig. 1. — Two examples of local a distribution observed in Hinode (SOT/SP) data. Back- 
ground is the continuum image. Red and blue contours represent positive and negative values 
of a respectively. The contour levels are ±1 x 10 _8 m _1 , ±5 x 10 _8 m -1 , and ±10 x 10 _8 m _1 
. The values of vertical current and a along the arc shown in the penumbra of the image in 
the left panel, is plotted in Figure 2(a) and those along the straight line in the umbra are 
plotted in Figure 2(b). 
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N0AA1 0933:05Jan07:0708hrs 

: (a) mean current= -0.896350 ± 41.0054 

= mean alpha= -0.887750 ± 50.9784 
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Fig. 2. — Plots of vertical current and a values along (a) the arc and (b) the straight line 
shown in the penumbra and umbra of AR NOAA 10933 (left panel of Figure 1) respectively. 
Black and red colors represent the current and a values respectively. The mean values of 
both the vertical current and the a values with their la standard deviations in the arc and 
the line are printed on the plots in their respective colors. 
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Fig. 3. — Cartoon diagrams of circular unipolar spots with positive and negative chirality 
are shown with the directions of observed transverse field (B obs ) and potential transverse 
field (B pot ). Solid and dashed lines represent observed and potential fields respectively. In 
this cartoon, the B f, s is tangential to the solid curved lines, while B pot is tangential to the 
dashed curved lines which have lesser curvature than the solid lines. First two cases bear 
positive chirality and later two negative chirality. Plus and minus sign in the central circular 
region represents the positive and negative polarity respectively. For details see the text. 



